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^ Abstract 
(N 

The Qth-power algorithm provides a useful canonical P-module presentation for 
the integral closures of certain integral extensions of P := F[x n , . . . , x{\, a polyno- 
mial ring over the finite field F := Z q of q elements. Here it is shown how to use this 
for several small primes q to reconstruct similar integral closures over the rationals, 
F := Q, using the Chinese remainder theorem to piece together presentations in dif- 
ferent positive characteristics, and the extended Euclidean algorithm to reconstruct 
small rational fractions to lift these to presentations over Q. 



C3 



£> . Key words: integral closure, normalization 

o 
o 

£T) ■ 1 Introduction 



The Qth-power algorithm, [6] , [5] , [I] , is designed to compute integral closures of 
integral extensions of multivariate polynomial rings in positive characteristic 
by exploiting the linearity of the Frobenius map a t— > a q in characteristic q. 
Over the rationals there is no such map to exploit. 

Since the integral closure over Q (the rationals) should specialize to the inte- 
gral closure over Z ? (the finite field of q elements) for almost all primes q, it 
makes perfect sense to compute the integral closure of the mod q image for one 
or more primes, use the Chinese remainder theorem to reconcile the results 
modulo the product of those primes, and then use the extended Euclidean 
algorithm to lift those results to results in characteristic 0. It will be proven 
that the presentation over Q is isomorphic to a subring of the presentation of 
its image over Z q in general (at least when everything involved makes sense 
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mod q), with equality for most q; and that the reconciled version of the pre- 
sentations over Z q for various q is lifted back to something over Q which has 
fractions with the same set of leading monomials as those of the fractions in 
its mod q images. So if the lifted version has an isomorphic image of the orig- 
inal ring inside it (which will not happen unless the product of the distinct 
primes used is large enough) , it necessarily must be the integral closure of that 
original ring. 



The fact that the extended Euclidean algorithm gives essentially inverse results 
of the mod q map when q is sufficiently large, gives, as a corollary, the integral 
closure for large primes q can be gotten from those for several smaller primes, 
a useful result here, in that the Qth-power algorithm should be, by its nature, 
expected to perform significantly better for smaller primes q. 



It should be pointed out that in both characteristic and characteristic q > 
the main advantage of the Qth-power algorithm is that it takes highly struc- 
tured input and produces highly structured output, namely a strict affine 
P-algebra presentation for the integral closure with an induced (as opposed 
to default) monomial ordering based on the weighted monomial ordering on 
the input. This allows for a fairly simple determination of the existence of a 
better Noether normalization than the given one, giving a smaller presenta- 
tion with the same type of structure and information. But in characteristic 0, 
it gives a presentation with relatively small rational coefficients rather than a 
presentation over the integers with overly large integer coefficients, and one 
that specializes mod q to that of the image mod q for all large q (and a subset 
of the integral closure for all smaller q for which the image makes sense). 



Section 2 contains notation to describe the main algorithm, the algorithm 
itself, and an outline of what will be proved to justify it. The technical def- 
initions and other details are postponed to section 3. Section 4 deals with 
the computation of canonical conductor elements based on the Jacobian. Sec- 
tion 5 is a discussion of the application of the Chinese remainder theorem 
and the extended Euclidean algorithm in this context. And section 6 contains 
the theory and proofs neeeded to justify the algorithm. There are numerous 
small examples throughout to help with the concepts and notation. But there 
is a larger example relegated to the Appendix useful as well for comparison 
to other implementations of integral closure computations. There are several 
other such on the author's website. Code for Qth-power algorithm has been 
available for a while in Magma and more recently in Macaulay2 as well. 
The latter has documentation containing more examples. 
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2 Overview of the algorithm 



Let P^ := Q[x n , . . . , Xi] be the polynomial ring over the rationals in the 
(independent) variables x n , . . . ,xi. Let := P^[y]/(f(y)) for some monic 
polynomial f(T) G P[T], be an integral extension. Suppose also that it is an 
affine domain of type I with a weight function defining a weight-over-grevlex 
monomial ordering (as in [6]), though more general hypotheses may suffice. 

Let Q(S W ) denote the field effractions of S {0 \ Let A< > G P (0) be the canonical 
monic conductor element computed from Jacobian(I?( )) as described below, 
so that the integral closure C(S (0) , Q(S (0) )) of in Q(S^) satisfies 

S ( ) Q c(S<®,Q(S®)) C -L^ ) C Q(S<®), 

and is known to be the union of all rings lying between and -^jS^°\ 
[Since is assumed to be an integral extension of p(°\ it probably makes 
more sense to think of it as C(P^°\ Q(S^)), in that P^ is a minimal subring 
over which this has a finite module structure. And once a fixed conductor 
element A^ G P^ has been chosen, it makes more sense to use the notation 
C(P(°\ -£jojS^) (even though ^oy^ '' is n °t a ring) to emphasize that the 
elements are from -^jS^- ' and are integral over p(°\ So that is the notation 
that will be used here.] 

The objective here is to find a canonical ordered set of monic polynomials 
(that is, with leading coefficient 1 relative to the monomial ordering being 
used) g^, . . . ,gi \go°' ) := so that the corresponding fractions gf^/A^, 
< j < Jo form a P^-module generating set for C(P^°\ -^jS^), then use 

y^jg , • • • , y^ as new variable names for the fractions g^ / A^°\ . . . , gf^/A^ to 
define a ring 

P (0) := Q^,...,^; £„,..., xi] 

with grevlex-over-weight monomial ordering induced by the weight function 
of (as in [6] but described below as well), and then compute the monic 

polynomials b£\ 1 < k < K , defining the minimal, reduced Grobner basis 

of the ideal 1^ of induced relations. 

Then^ (0) :=E (0) /7 (0) is a strict affine P^-algebra presentation of the integral 
closure C{P^\ -^S^). While is an explicit subring of S { °\ need 

not be; so let : — > be the inclusion map (the identity on the Xi 
identified in each copy of P throughout, but not on the y which is mapped to 
combinations of the y~j and the Xi) so that 

^(o)(5-(o)) g 5 (0) C -L^)^)) c Q(ipP\sV>)). 
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Now consider what needs to happen for there to be an image of all this over 
Z q for some prime q, gotten by applying the mod q map, \i q . 



It is easy enough to define P^ := Z q [x n , . . . ,xi] by identifying the variables 
Xi, and similarly to define := Z q [y;x n , . . . ,x±] by further identifying the 
variable y. 

If q doesn't divide any denominator j3 of any rational fraction a/ (3, gcd(a, (3) = 
1 of any fcf, 1 < k < K, then := := ix q {bf) G R (q) : 1 < k < K) 
makes sense, and is still a minimal, reduced Grobner basis for the ideal 1^ 
of R^ that it generates, though the quotient ring 

S ( q ) ;= #(?)//(<?) nee( i 

no 

longer be even a reduced ring, let alone an afline domain of any sort. [Whether 
it is an integral extension of P^ can depend on how one views such things as 
Z q[y i'i x i]/ (Vi) m which xi doesn't appear in the defining relation.] 

And if q doesn't divide LC(d ( V), for = LC(d^)A ( ^ the conductor as 
computed from the Jacobian of R>( > over Z, then A(«) := fi q (A^) is the 
monic canonical conductor element that would have been computed using the 
Jacobian over Z q . 

The Qth-power algorithm is meant to work in positive characteristic to pro- 
duce a strict afline P^-algebra presentation with P^-module generating set 
of fractions with monic numerators g^ / A^ q \ . . . , g[ 9 ^ / A^ q \ and gj^ / A^ := 1, 

given the variable names y^, . . . ,y^ to define 

S (9) :=F ff |^>,...,#;x B ,...,x 1 ] 

(having the grevlex-over- weight monomial ordering induced by that on R^) 
with monic polynomials b k , 1 < k < K q forming a minimal, reduced Grobner 
basis B for the ideal I of induced relations, defining the presentation 
S (q) :=R {q) /f q) ■ 
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The steps in the proposed characteristic algorithm based on this are then 
simple to understand: 



Algorithm 1 (1) Start with the finite ordered set of {independent) variables 
(x n , . . . ,Xi) defining the Noether normalization P^ in characteristic 0, 
the ( dependent) variable name y used to define the ring R^°\ and the finite 
ordered set of monic relations (b±, . . . ,bx) forming a minimal, reduced 
Grobner basis for the ideal of relations 1^ for a presentation of the input 
quotient ring = R^/I®. 

(2) Compute a canonical conductor element A^ G P^ for from the 
Jacobian . 

(3) For successive primes, qi, test that the (mod q{) map, fi qi , is defined (that 
is that qi doesn't divide (3 for any rational coefficient a/ (3, gcd(a, (3) = 1, 
in any of the basis relations bk or of A^) and that is still an integral 
extension of the image 

(4) Compute a canonical conductor element A^ G P^ for S^ qi \ skipping 
qi if it is one of the (finite number of) primes for which A^ ^ fi qi (A^). 

(5) Use the Qth-power algorithm in characteristic qi (as a black box for the 
purposes of this paper) to compute a canonical ordered set of (numerator) 
polynomials (g^), < j < J qi (with the common denominator polyno- 
mial being g^ := A^)) for the fractions forming a P^ -module generat- 
ing set for a strict affine P^ -algebra presentation of the integral closure 

C ( p{qi) i^h siqi) ) with (bk l) ) , 1 < k < K qi , the minimal, reduced 

Grobner basis for the ideal of relations relative to the induced grevlex- 
over-weight monomial ordering. 

(6) If (q{), I E L is a sequence of distinct primes for which presentations 
have been computed, and 

• Jl '■= Jqi is independent of I G L; 

• LM(g^) is independent of I G L; 

• Kl '■— K qi is independent of I G L; 

• LM(b^) is independent of I G L; 

then use the Chinese remainder theorem on the canonical ordered sets 
(?j 9 ^\ eL > 1 — J — '^ L ' an< ^ a ^ so on ^ e se 1 uences i^k^j , 1 < k < Kl, 
to get similar canonical ordered sets (gf 1 ^), 1 < j < Jl, and sequences 
T>^ L ), 1 < k < K L , for N L :— I\{qi '■ I G L} monic with the same sets 



of leading monomials. 

(7) Then use the extended Euclidean algorithm to lift the coefficients c(mod N L ) 
to small fractions a/ (3 G Q with a 2 + (3 2 minimal, to get a canonical or- 
dered set of polynomials (gf ,Nh ^), 1 < j < Jl, and sequence (b^^^), 

1 < k < Kl, over the rationals with the same sets of leading monomials 
to describe a possible integral closure in characteristic 0. 

(8) Stop when {bf NL) ), 1 < k < K L , is a minimal, reduced Grobner basis 



5 



for the ideal I ' L it generates and the image ^ Nl \I^) of the original 
ideal 1^ is contained in I^°' Nl ^ ; both necessary conditions. 

To explain why this works, the important steps are to show that whenever 
everything involved makes sense 

(i) 

(2) S^°' Nl ^ i s isomorphic to if 1^ is isomorphic to a subideal of I^°' Nl ^ 
and B^°' Nl ^ := (frj, '^), 1 < k < K L , is still a minimal, reduced Grobner 
basis for I^°' Nl \ 

These are proven as Lemma 9 and Theorem 15 below. It should be noted that 
(regardless of characteristic) any P-module between S and -^S has a canonical 
ordered set of polynomials (gj)j (as defined in the next section) with {LM(gj)j) 
a measure of the size of the P-module. That is, were the de Jong algorithm 
implemented relative to a fixed Noether normalization P and a canonical 
conductor element A e P, then the sequence of nested rings produced would 
have a sequence of canonical ordered sets of polynomials (gj)j with (LM(gj)j) 
nested and getting larger. The reverse is true of the Qth-power algorithm 
in that the P-modules produced have canonical ordered sets of polynomials 
with (LM(gj)j) nested and getting smaller. Both approaches meet in 
the middle with a P-module that is a ring that must be the integral closure 
sought. 

Also, Jacobian(p(°)), over the integers, Z, is used below to define canonical 
conductor elements A^ ^ e P^ and e P^ for all primes q at the same 
time. This is discussed in its own section. The use of the Chinese remain- 
der theorem and the extended Euclidean algorithm, while discussed below as 
applied in this context, are assumed to be elementary. The proof that the 
Qth-power algorithm works in positive characteristic was dealt with in the 
author's previous papers cited in the introduction, though certain parts of it 
are discussed below. [It should be noted however, that it would not take too 
much work to put other implementations of other integral closure algorithms 
in a form that would also work here, though at present, few if any give a sim- 
ilar canonical result in characteristic that directly specializes to the result 
they give in positive characteristic. It makes mathematical sense to rewrite 
them to reflect this connection between integral closures in characteristic 
and positive characteristic q.} 
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3 Definitions and other details 



The following material describes the structure that is used to describe inte- 
gral closures of integral extensions of a given Noether normalization P, and 
explaining the mindset of the Qth-power algorithm approach to same. The 
idea is to have an integral extension S :— R/I of P := F[x n , . . . ,xx), with 
R := F[y; x n , . . . , xi], I :— (f(y)) the ideal of relations, and P a Noether nor- 
malization of S, and with S having a weight function induced on it by P, [B] . 
Then its integral closure C(P,S) should have a presentation S := R/I with 
an induced weight function. That is, for R := F[y~j, ■ ■ ■ ,y 1 ,x n , . . . , xi] with 
y a name for the (non-trivial) P-module generator gj/go having wt(y~j) : = 
wt(gj) — wt(g ) as its induced weight. 

Moreover a presentation S of the integral closure C(S,Q(S)) should have a 
nice structure as an affine P-algebra. 

Definition 1 A strict affine P-algebra presentation S := R/I with R := 
F[y, x], and I the ideal of induced relations is one with a minimal, reduced 
Grobner basis B for I consisting of P -quadratic relations of the form y-y^ — 
J2k c i,j,kVkJ c i,j,k £ P, describing the P-algebra multiplication with possibly 
some monic P -linear relations of the form J2k a kVk! a k €= P, if the P-module 
generators, y k , are not independent over P, 

This is ensured by the grevlex-over-weight monomial ordering, in that all 
products of total degree 2 in the y's are reduced to P-linear combinations of 
the P-module generators with total degree less than two, and all P-syzygies 
only involve monomials of total degree less than 2 in them. 

Definition 2 For P := F[x n , . . . , Xi] a polynomial ring, R := F[y; x n , . . . , x±], 
I an ideal of R such that the quotient ring S := R/I is an integral extension of 
P, an ordered set (gj G R : < j < J) of polynomials is said to be canonical 
for some submodule 4T C ~S iff 

(1) each gj is monic (has leading coefficient 1 relative to the monomial or- 
dering being used); 

(2) go = A G P is a conductor element for T; 

(3) (gj/ 'A : < j < J) is a P-module generating set for ; 

(4) (gj '■ < j < J) is interreduced, meaning that no monomial of any gj l 
is of the form x—LM(gj 2 ) for any ]2 ^ ji and x— G P. 

Making sure that one produces a canonical strict affine P-algebra presentation, 
independent of characteristic, is crucial in being able to reconstruct a canonical 
presentation in characteristic from canonical one in positive characteristic. 

Moreover the induced grevlex-over-weight monomial ordering is based on a 
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weight function on the input extended to the output: 

Definition 3 A weight function wt : R\I — > N n (n the number of indepen- 
dent variables) is a function satisfying: 

(1) wt(c) = for all c E F\{0}; 

(2) wt(gh) = wt(g) + wt{h) for all g, h; 

(3) ifwt(g) = wt{h), then for some unique c G F either wt(g-ch) -<i ex wt(g) 
or g — ch E I . 

If R — P — F[x n , . . . , Xi], the columns of any non-singular matrix Mp defining 
a (global) monomial ordering on P also define a weight function on (the non- 
zero elements of) P. 

Example 4 F[y,x]/(y 3 — x 5 — yx) has a weight function with wt(y) := 5 
and wt(x) := 3 (corresponding to pole orders of these as rational functions at 
the only point where either has a pole). The example F[y, x]/{y 3 — x 3 y — x) 
does not. [An attempt at weights would have wt(y 3 ) = wt(x 3 y), so 3wt(y) = 
wt(y) + 3wt(x) > wt(x) > with no field element c with wt(y 2 — cx 3 ) a 
smaller weight because y 2 — x 3 = x/y would have weight less than 0.] The 
example F[x, y, z]/ (x 6 + x 3 z — y 2 z) has a weight function with wt(x) := (3, 2), 
wt(y) := (6,6), wt(z) := (6,0), an extension of the grevlex monomial order 
with wt(y) = (1, 1) and wt(z) = (1,0). 

Weight functions (relative to a given ideal I) have the important property that 
wt(g) = wt(NF(g, I)), as otherwise their difference (an element of I) would 
have a defined weight. This, in turn, implies that all standard monomials have 
different weights. Integral extensions with weight functions have at least this 
much more structure than those that don't. The integral extensions such as 
those considered here have a weight function. 

One can extend a weight function on S naturally to (the non-zero elements 
of) Q(S) (the field of fractions of S) by wtQ^s){g/h) := wts(g) —wts(h) if one 
allows values in Z n . But for g/h £ C(S, Q(S)) : wt(g/h) will not have negative 
entries and will represent an induced weight on g/h. 

There are weight-over-grevlex and grevlex-over- weight monomial orderings de- 
fined by the non-singular matrices gotten by replacing the top n or bottom n 
rows of a grevlex monomial ordering matrix by the weight matrix, respectively 
The former emphasizes the property that wt(LT(f)) = wt(LT(f — LT(f))), 
whereas the latter emphsizes the desired strict affine P-algebra presentation. 

The conductor element A and the numerators gj produced are all assumed 
to be monic. What will be computed are pairs of finite canonical ordered 
sets of polynomials and finite sequences of relations, with maps between such 
pairs. The induced weights are kept track of as well, given that they define 
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the monomial orderings involved. 



What is necessary to know about the Qth-power algorithm is that it treats 
the input ring S as a P-module with a natural induced monomial ordering, 
computes a conductor element A e P, starts with a dual module such as the 
default M := j^S and computes a nested sequence of P- modules 

M D • • • D M L = M L+l = C(P, ±S) 
by the simple definition 

M m :={£eM i: NormalFarm^,!) 

Necessarily each Mj (and hence the integral closure itself) is a P-module with 
a natural induced monomial ordering. Moreover it naturally produces a strict 
afline P-algebra presentation R/I relative to a canonical ordered set (gj)j 
of polynomials. [This approach works theoretically for any characteristic and 
any integer power at least 2, but is only linear when q is (a power of) the 
characteristic! 



Consider the reasonably generic example: 



Example 5 Let := Q[x), Pi ) := Q[y;x], and J(°> := (B^), with 

p(°) : = {( y 2 - 3/Ay - 15/17x) 3 - 9yx 4 (y 2 - 3/% - 15/17s) - 27x u }, 

with weight function (11,6) (meaning that the monomials y 6 and x 11 are the 
only two in the defining relation of largest weight, 66) defining an integral 
extension S<® = P (0) // (0) - [This was originally constructed to have at least 
one non-trivial integral element, namely (y 2 — 3 /Ay — 15/17x)/x 2 and several 
moderate- sized rational coefficients to be reconstructed, in order to test this 
extension of the Qth-power algorithm to characteristic 0.] 

This is worked out in detail in the Appendix, not only by the methods be- 
ing described here and implemented in the author's QthPower package in 
Macaulay2 [TO]/ , but also using the other existing applicable implementations 
of integral closure and/or normalization algorithms, normal in Singular, in- 
tegralClosure in Macaulay2, and both Normalisation and IntegralClosure in 
Magma. 
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The output from the qthlntegralClosure function in the QthPower package 
consists of: 

(1) an ordered set of "numerators" for a P(°) -module generating set {here a 
basis) : 

(0) / 4 2 3 ^ 3 15 4 

g \>/x = y x - -yx - —x , 

(0) / 4 5 

g\ ' jx = yx , 

(o) / 4 4 3 o 30 no 9 2 45 2 225 o 

Oo ' x =yx yx -y x H y x H -yx H x, 

y 3 / y 2 y 1? y lg y 34 y 2gg > 

(o) ; 4 s s 15 4 9 o 45 4 
a\ x = y x yx yx x , 

(0) . 4 . 9 4 30 o 27 o 45 2 225 2 27 2 135 675 2 

= V — "V V x-\ y H y x-\ yx y yx x 

y A y \T 16 y 17 y 289 y 64 y 136 y 1156 

with g^ /x 4 = 5^ = x 5 e P^ also being the {reduced) conductor element 
and common denominator of the fractions, {not the actual = x 9 , 
accessible only through the qthConductor function); 

(2) a minimal, reduced Grobner basis of induced relations; 

( 3) an induced weight matrix {here a weight matrix) (25, 21, 20, 11, 10, 6) defin- 
ing a grevlex- over- weight monomial ordering. 

A typical induced relation in the Grobner basis, such as 

V4V2 ~ 3/2y 4 - y z x\ - 15/Uy^ + 9/8y 1; 

{produced by finding the normal form ofy 4 y 2 relative to the input ideal 
corresponds to some P-algebra multiplication rule, in this case 



/21 * fn = 3/2/ 21 + f 6 f 20 + (15/17/ 6 - 9/8) / 
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{using notation f 25 : = y 5 , /21 := y A , /20 := V 3 > fu '■= V2> fw '■= Vi, h ■= x i to 
reflect the weights). Note that the weight of the left side, /21/n and the weight 
of its NormalForm, the right-hand side, are both necessarily the same {in this 
case, 32). 
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4 Computing a canonical conductor element 



Standard methods to compute a conductor element A G S (meaning an el- 
ement for which C(S,Q(S)) C —5) use determinants of n x n minors of 
a Jacobian. This can be easily done by column-reducing the Jacobian ma- 
trix Jacobian(B) of B; and this computation can be done over R instead 
of S by appending columns one for each basis element of I and each row 
of Jacobian(B). It is then possible to consider those entries Cij G P of the 
column- reduced form, C, for which C& j = for k > i and I < j. An appropri- 
ate monomial ordering must be chosen relative to which this is done, so that 
the elements Cjj G P considered will correspond to diagonal entries innxfi 
minors whose determinants necessarily produce conductor elements, greatest 
common divisors of those in the same row can be used, and a scaled product 
of those gcds over all rows can be used to give a canonical conductor element 
A G P. For this purpose, any block ordering treating the dependent variables 
any way but using the given monomial ordering described by Mp on the lowest 
block consisting of the (independent) variables in P, will suffice. [Note that 
when computing A^ ^ G P^ over Q, it is possible to do this over the integers, 
Z, instead (if denominators are cleared first) in order to see in one computa- 
tion for what (finite set of) primes q it might be that ^ A^(mod q), by 
seeing what primes occur anywhere in the column reduction C] 

The method qthConductor exported from the author's QthPower package, [TO] , 
in Macaulay2 can be used to compute such a canonical conductor element, 
by letting Macaulay2 do the column-reduction, then using a simple loop 
to compute the product of the gcds described. This computation is not a 
point of this paper, other than to insure that there is a canonical conductor 
element that can be computed, that it is an element of the given Noether 
normalization, P, and that the computation in positive characteristic mirrors 
the computation in characteristic 0. 

Consider the following instructive example, meant originally to test minimality 
and form of presentation, but, as a byproduct, was used to catch bugs in 
various implementations as well. 

Example 6 Let p(°) := Q[x], / (0) (T) := T 8 - T 2 x 3 + 2Tx 6 - x 9 G P^[T], 
R(°) := Q[y, a;] with a monomial ordering based on W$ '■= (9,8), and := 
R(°)/(f(°)( y )). The (extended) Jacobian matrix (df (0) (y)/dy, df i0 \y)/dx, f {0) (y)) 
column-reduces over Z with lex, x -< y, monomial ordering to 

(120a; 24 , 15a; 26 , 5x 27 , a; 41 + 6a; 26 , 

6yx 5 - x 38 - 24x 23 - 6a; 8 , 2yx 6 - 2x 39 - 48x 24 - 2a; 9 , yx 26 + 2x 29 , 
3y 2 x 2 - 2x 38 - 48a; 23 - 3x 8 , y 2 x 6 - 2x 27 - a; 12 , 8y 7 - 2yx 3 + 2x 6 , 
y 7 x 2 + yx 20 + 2yx 5 + x 38 - 10x 23 - 2x 8 , y 8 - y 2 x 3 + 2x 39 + 48a; 24 + x 9 ). 
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From this it is easy to extract a conductor element A(°) := x 2A G P^ . It is 
also easy to extract A® := x 2A E P (q) for q^2,3, 5, A (2) := x 26 , A® := x 27 , 
and A^ := x 26 (x 3 + l) 5 for similar problems with F = Z q , q a prime. 

There are rational functions (f := 1, / 4 , / 5 , /g, fw, fu, /15, /19) (with the sub- 
scripts corresponding to the weights) forming a P^ -module basis for the in- 
tegral closure, . Then the presentation of can be gotten by reading 
modulo q for all primes q 7^ 3,5. Curiously, the smallest conductor el- 
ement that could be used is 5^ = 5^ = x 13 for all primes except 5^ = 
x 13 (x 3 + l) 2 . It is tempting to conjecture that A^ = A^(mod q) implies that 
g(q) _ gW^jnQfi qy jf ^ c \ ea rly not true that 5^ = 5^ (mod q) implies that 
g(q) _ g^^mod from q = 3 in this example; and it is clearly not true that 
S (q) = S (0) (mod q) implies that A^ = A^(mod q) from q = 2. 

Since it is computationally easy to avoid all the (finitely many) primes q 
for which A^ 7^ A^(mod q) (necessarily divisors of some coefficient in the 
computation over Z), it is possible to simplify subsequent computations by so 
doing. 
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The exportable QthPower code in MACAULAY2 for this is [TO] : 

Algorithm 2 qthConductor = method (TypicalValue => RingElement) ; 
qthConductor ( Ideal ,ZZ) := (I,depno) -> ( 
R := ring I; 

RP := (coef f icientRing R) [gens R,MonomialOrder=> 

{Position=>Up,-[depno,#gens R-depno}}] ; 
IP := sub(I.RP); 

GP := gens gb (transpose jacobian IP I 

matrix-[{gens IP}}-**identity(RP~ (numColumns (jacobian IP)))); 

depvars : =take (gens RP , depno) ; 

indvars:=take(gens RP , depno-#gens RP) ; 

qthconductor := 1; 

rowconductor := 0; 

j := numColumns(GP)-l; 

i := numRows(GP)-l; 

while i >= and j >= do( 

while i>=0 and j>=0 and (GP_(i,j) == 

or (logpoly(GP_(i , j) , depvars , indvars) )#1 != 1) do( 

j = j-i; 

); 

rowconductor = 0; 

while i >= and j >= and GP_(i,j) != 
and (logpoly(GP_(i, j) , depvars, indvars))#l == 1 do( 
rowconductor=gcd (rowconductor ,GP_(i , j) ) ; 

j = j-i; 

); 

if rowconductor != then qthconductor = qthconductor*rowconductor ; 
i = i-1; 

); 

s : =sub (qthconductor , R) ; 
s/leadCoef f icient (s) 
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5 Chinese remainder theorem and extended Euclidean algorithm 

As stated above, the approach given in this paper is, in some sense, an elemen- 
tary one, in that it exploits commonly known information from the Chinese 
remainder theorem, and intermediate information given by the extended Eu- 
clidean algorithm. Even a good reference such as [2] doesn't necessarily use the 
extended Euclidean algorithm in this exact way. And, as a warning, this ap- 
proach is extremely tricky in the sense that the maps are not homomorphisms 
of the whole rings involved, but do extend naturally to ring homomorphisms 
when defined correctly on the finite ordered sets of objects used to describe 
those rings. That is, a presentation of the integral closure of an afline P-algebra 
can be described in terms of a finite ordered set of fractions and a finite or- 
dered set of relations. It is then possible to construct those two finite ordered 
sets, define mappings, and then extend those naturally to definitions of ring 
homomorphisms. 

Note especially that in implementing this approach, care must be taken to 
assure that the integral closure algorithm produces the same canonical result 
for each good prime q. That is, for most primes q, the integral closure over Z q 
should look exactly like that of the integral closure over Q, but with coefficients 
reduced mod q. 

For each presentation of S := R/I and presentation of its integral closure 
S := R/I, there is a map ip : R — >■ R, necessarily with if>(I) C J, so that ip 
can be viewed as an inclusion map ip : S — > S. 

The extended Euclidean algorithm can be used to move between % e Q, with 
gcd(at, j3) = 1, (3 > and representatives c 6 Zjv. 

Definition 7 The rational fraction reconstruction map (see, for instance, JMf) 
is 

E N (c) := — , c/3 = aimod N), a 2 + (3 2 minimum , (3 > minimum. 
The mod N map is 



These are almost inverse operations in the sense that for —N/2 < c < N/2, 
fi N o E N (c) = c; while, for a 2 + (3 2 < N, E N o fi N (%) = |. 



Both maps naturally extend to polynomials, by applying them to coefficients 
and mapping variables to corresponding variables; so we shall abuse notation 
and use the same function names when applying them to polynomials. 




c, c/3 = a(mod N), |c| minimum . 
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Definition 8 Similarly the Chinese remainder map standardly takes ordered 
sets of remainders (a/); e L an d ordered sets of respective moduli (qi)i^L, and 
produces a(mod N L ) for N L := l\{Qi '■ I £ L}) such that a = aiiynod qi) for 
all I G L when the moduli are all relatively prime, as they will necessarily be 
here when the qi are distinct primes. 

We shall call this map CRT regardless of the number of inputs, and regardless 
of whether we are applying it to integers or extending it to polynomials by 
applying it to the coefficients. 

Lemma 9 fi q (C(P (0) , -^ojS^)) C C(P^\ -^S^) for all primes q for which 
fi q (C(P(°), ^ {0) )) makes sense. 

PROOF. For each fraction gf^ /A^ in the desired integral closure over Q, let 
ff\T) G P (0) [T] be a monic polynomial satisfied by it. If both fi q (gf) / n q (A<°>) 

and ff\T) := fi q } [T]) G P^[T] are defined (meaning the the prime q 
doesn't divide the denominator (3 of any rational fraction a/f3, gcd(a,(3) = 1 
occurring in either gf or /j 0) [T]), then fx q (^ 0) ) / fi q (A( >) G ^jS^ satisfies 

the monic polynomial f- q) (T) G P {q) [T}. 
Definition 10 A prime q is a good prime iff 

^(^ P(0) '^ (0) ))=^ P(9) '^ (9) )- 

Corollary 11 // (gf : < j < J) is canonical for C(P {0 \ ^jjS^) and the 

mod q map makes sense, then (gf ] : < j < J) is canonical for C \P^\ ^jS^) 
if q is a good prime (and only for some subring if it is not a good prime). 

Clearly if q divides any denominator of any rational coefficient a//3 of any 
term of any bf\ it is not good. If q divides any numerator of any rational 
coefficient a//3 of any term of any bf \ it is may not be good, especially if the 
extension mod q is no longer really an extension. And if A^ ^ fi q (A^ ^, q 
may not be good. So computationally one can try to avoid such primes that 
are not good or may not be good (since these form a finite predictable set of 
primes) . 

Example 12 Consider the example with £? (0) := {yf + 13/22(x^ + x\ + x\)}, 

for which we should expect y\ = y\jx\ and B^ = {y\ + 13/22 (xf + x\ + xi)}. 
The primes q = 2, 11 are clearly bad since the mod q map, /j, q , makes no 
sense; but q = 13 is also bad in the sense that = {yf} really doesn't 

define an integral extension of P^ = Z 13 [a;i]. Column-reducing the Jacobian 
(22yf + 13 (x I +xl+xl ), 44yi, 13(9rrf + 7x\+hx\)) overZ with lex ordering gives 
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(780:4 156xf , 52x? + 260xf , 13x? + 39x? + 65x?, 4% b 22y 2 + 26x[ + 52xf ) . So 
£/ie canonical conductor elements are A^ = A^ = xf /or a// oi/ier g except 
^(3) _ ^.6 _ x 4 an( ^ ^(5) - ^ g _ 5 happens to be a good prime in this 

example, but for q = 3, y 1 = yi/(x\(x\ — I)), meaning there is a larger than 
expected integral closure S . Avoiding the primes 2,3,5,11,13 {whether or 
not they are not good), using 7,17,19 (which should be good) is enough to 

reconstruct S {0) , since 7 • 17 • 19 = 2261 > 22 2 + 13 2 = 653. 

The Euclidean algorithm, applied to Nl '■= r_i and any r > 0, produces 
sequences (rj) and (Qi) such that rj_ 2 = QiU-i + r 'i with < r\ < rj_i, and 
r n = 0. Part of the extended Euclidean algorithm produces a sequence (ui) 
with U-i := 0, u := 1, and U; L := + Wi_2- Then for each i, (— l)Vj/iij = 

ro(mod Nl). Of these there is necessarily some % > with r 2 + -u 2 minimum, 
choosing % minimum as well if this is not unique. 

Now define the composite map 

V (W := E Nl o CRT o :leL})o : I e L}) 

for ip( q ri the corresponding inclusion map from to if . 

Suppose the variables in the integral closure presentation correspond 
to the fractions gf^ / A^ for gf \ e i?^, and the Grobner basis elements 
of the ideal of induced relations I are denoted by b k . Let the variable iff 1 
correspond to gf/A^ q \ for gf := /x 9 an d := (At )). If q is a 

good prime, then these should be variables and (a Grobner basis of) relations 
for the integral closure S^ 9 \ 

Since the object here is to go in the reverse direction by reconciling various 
presentations, S^ q \ and reconstructing the presentation from them, us- 
ing the Chinese remainder map and the extended Euclidean algorithm map, 
consider the candidates for S^°\ namely 5 ,< -°' 7Vi - ) with polynomial ring ^°' ArL - ) 
having variables yf' NL ^ corresponding to gf ,Nl ^/A^ for 

gf NL) :=E NL oCRTo\{{g^: I £ L} 
and ideal / < -°' 7Vi - ) generated by the finite ordered set of images 
t NL) ■= E Nl o CRT o n{6? } : I EL}. 
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6 Theory 



Lemma 13 If q — Ni is a good prime larger than a 2 + f3 2 for any coefficient 
a//3 e Q needed to be reconstructed to produce the presentation / then 
j^ii) fojfe ^ Q p resen i a n on ^ [And the canonical polynomial set (g^ : < 
j < J) necessarily lifts to a canonical polynomial set (gf' q ^ : < j < J).] 



PROOF. If gf lifts to gf (including g$> = lifting to ^ 0) = A (0) ), 

and the relations bf lift to bf\ then S (q) lifts to S W . But if q > a 2 + (3 2 
then c = a/ (3 (mod q) lifts to a/j3 using the extended Euclidean algorithm as 
described above. 



Corollary 14 If (qi)ieL is a set of distinct good primes and Nl '■— Y\{qi '■ I € 
L} is larger than a 2 + (3 2 for any rational coefficient needed to be reconstructed 

to produce the presentation /I^°\ and R^ / 1^ are compatible in the sense 

that LM(g^) is independent of I and LM(b^) is independent of I, then 

(R^/I^)i e L lifts to this presentation. [And, again, the canonical polynomial 
set (gf 1 ^ : < j < J) necessarily lifts to a canonical polynomial set (gf' NL ^ : 
0<J <■!)■] 



PROOF. Use the Chinese remainder theorem to reconcile these individual 
presentations, and lift the resulting ordered sets (g^^), 1 < j ' < Jl, and 

(b[ NL) ), 1 < k < K L , to ordered sets (gf' NL) ), < j < J L , and C^ '^^, 
1 < k < Kl, and proceed as in the previous proposition. 



Since is not known ahead of time, it is not clear how big N L must be 
to apply the proposition or corollary above. It is therefore better to have a 
theorem independent of this knowledge. So the following is a way of knowing 
that Nl is sufficiently large without knowing just how large sufficiently large 
is. 

For S^ '^ to be a presentation of the integral closure of , it necessarily must 
be a ring containing and also contained in -^jS^, being (isomorphic 
to) the union of all such. 

Theorem 15 IfB ( °' N) is a Grdbner basis for • 7 (M0 , and ^ N \l^) C/ (0,JV) , 
thenS ( °' N) =^(S {0) ). 
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PROOF. If B { °' N) is a Grdbner basis for I { °' N) , then the quotient ring S ( °' N) := 

j£0,N) ^j(0,N) . g a gtrict affine p(0). algebra> 

If ^W(j(°))c7 (M0 , then 

^) (S (»))cs M clf)(S( '). 

But V (iV) (^ (0) ) is the union of all such rings, so S { °' N) C ^)(S (0) ). 

If S (0,A0 ^ ^^)(5 (0) ), consider the monic conductor element A^ G P (0) 
mapping to A^ e P^ for all primes g not identified as bad primes by the 
Jacobian computation above. Were the integral closure of tp( N \S^) computed 

as the integral closure of S^ '^, the conductor element A(° ,Ar) e P^ computed 
would necessarily be a divisor of A^ ^ since S {0 ' N) c ^(S^.SoA^) would 
be monic with mod q image n q (A^°' N ^) e for any ?|iV. But jj q (S' y0 ' NS) ) = 
for these g. Since S'*' 9 '' is integrally closed, /j, q (A^ ' N ^) = 1. Hence A^ '^, 

being a monic element of P(°\ would be 1 as well; meaning that S {0 ' N) would 
be integrally closed. 

Example 16 Let 

S (q) := F q [y 7 ; x 3 ]/(y 3 7 + x\ + 8y 7 x 3 ), 
(with Fq : = Q allowed). Then 

S (0) = Q[z n ,y 7 ; x 3 ]/ (z 2 n + y 7 xl + 8z n , z u y 7 + x\ + 8y 7 , y 7 -z n x 3 ), 

with the subscripts defining an induced weight function and a correspond- 
ing induced weight- over- grevlex monomial ordering. For any prime g ^ 2, 

fi q (^S^^j = S^ q \ So, for any N a product of distinct odd primes, = 

if N > 8 2 + l 2 For a smaller N such as N = 55, 8(mod 55) lifts to possibly 
the wrong fraction, here 1/7 instead of 8/1, and subsequently 

V 7 + X l + §2/7^3 = Vi{y 7 - zux 3 ) + x 3 (z u y 7 + x% + \y 7 ) + ^-y 7 x 3 £ J (0 ' 55) . 



For q = 2, 

S {2) = F 2 [w 1 ;x 3 ]/(w 3 1 +x 3 ), 
with y 7 = W\x\ and z\\ = w\x\. Clearly this is larger than expected, so S D 

s {0) . 
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Example 17 The generic example F q [yi] x\\/ (yf + a q y\X\ + b q x\) has y\/x\ 
in its integral closure; so its integral closure has a presentation as 

F q [y 2 , Vi, xi]/ (yl + a q y 2 + bq^xf, y 2 y 1 + a q y 1 + b q x\, y\ - y 2 x 1 ). 



If F = Q, and a := 1/3 and b := 8/7, then the image in characteristic 
q is not defined for q = 3,7, and is not an affine domain for q = 2. For 
F 5 = Z 5 , a$ = 2 and b$ = —1 would lift to a = 2 and b = — 1, giving a 
presentation of the wrong integral closure (one with the right form but these 
wrong coefficients). Using F n = Z n as well would give an = 4 and bn = —2 
reconciled to give a 55 = —18 and 6 55 = 9, and lifted to a = 1/3 and b = —1/6, 
again giving a presentation of the wrong integral closure. Using in addition 
F 13 = Z 13 would produce a± 3 = —3 and b 13 = 3 reconciled to give a 7 i 5 = —238 
and 6715 = —101, lifted to the correct a = 1/3 and b = 8/7. 

Example 18 Q[y, x]/(y 2 - 3/2x 3 + 2A/7x 2 - 96/49x) doesn't need N > 96 2 + 
49 2 to work, but only N > 8 2 + 7 2 , since the only things needed to be computed 
are go := x — 8/7, g\ := y, and ~b~i := y 2 — 3/2x (and the inclusion map image 
il>(y) :=y(x-8/7)). 

The details for this example are as follows: The primes q = 2,7 are bad because 
they divide denominators of fractions defining the problem. The image for q = 
3 is not even a reduced ring, so probably should be avoided as well. 5^ = 7x—8, 
and q = 7 is already to be avoided. 

Forq = 5, ^ 5) = A( 5 ) =x + l, g{ 5) =y, b {5) = y 2 + x, ip {5) (y) = y(x+l). This 

lifts to give g { °' 5) = x+1, g^' 5) = y, &[°' 5) = y 2 +x, ^(y) = y(x+l). Then the 
defining relation above reduces to —x(x + l) 2 — 3/2x 3 + 24/7x 2 — 96/49:r 7^ 0. 

For q = 11, g( ir > = = x + 2, g{ n) = y, = y 2 - Ax, ^ u) (y) = 

y(x + 2). This reconciles with the previous to get g^ = x — 9, g[ 55 ^ = y, 
5 (55) = y 2 + 26x, V (55) (l/) = V(x - 9). This lifts to give g^ 55) = x + 1/6, 
^(0,55) _ ^ ^(o,55) _ -2 _ 3/2x, ^( 55 ) = y(x + 1/6). Then the defining relation 
above reduces to 3/2x(x + 1/6) 2 - 3/2x 3 + 2A/7x 2 - 96/A9x ^ 0. 

For q = 13, ^ 13) = = x - 3, g{ 13) = y, b (13) = y 2 + 5x, ^ 13 \y) = 

y(x — 3). This reconciles with the previous to get g^ 1 ^ = x + 101, g^ 15 ^ = y, 
5 (715) = y 2 + 356x, V> (715) (l/) = V(x + 101). This lifts to give g ( °' n5) = x - 8/7, 
^(oji5) _ ^ ^(0,715) _ _ 2 _£^2x, -?/;( 715 ) = y(x — 8/7). Then the defining relation 
above reduces to 3/2x(x - 8/7) 2 - 3/2x 3 + 2A/7x 2 - 96/49x = 0. 

The presentation found (but not minimized) is then 

Q[y;x]/(y 2 -3/2x) 
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with inclusion map defined byipiy) = y(x — 8/7). [ The minimized presentation 
here would have been just the polynomial ring Q[y] with x = 2/3y 2 and y = 
y(2/3y 2 — 8/7) both unnecessary except for defining the inclusion] 
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7 Appendix 



It is envisioned that the code and the relevant examples relative to this paper 



on the website http://www.dms.auburn.edu/~leonada. will be updated as 



various packages change for the better. The code for the Qth-power algorithm 
in positive characteristic and the extra code to extend it to char for this 
paper are both written in Magma and in Macaulay2 and are available 
from the author. 

But, as mentioned in the Overview section above, the complete version of the 
example mentioned there, is done here by the various methods mentioned. 

P(°) := Q[x], R {0) := Q[y,x], 

: = {(y 2 - 3/Ay - 15/17x) 3 - 9yx 4 (y 2 - 3/Ay - 15/17x) - 27a; 11 } 



The Qth-power algorithm implementation produces fractions with numerators 
p_5:=x~5, 

p_4 : =y~2*x~3- (3/4) *y*x~3- (15/17) *x~4 , 
p_3:=y*x~5, 

p_2:=y~4*x- (3/2) *y~3*x- (30/17) *y~2*x~2+ (9/16) *y~2*x+ (45/34) *y*x~2 

+(225/289)*x"3, 
p_l:=y"3*x"3-( 15/17) *y* x ~4- (9/16) *y*x"3- (45/68) *x"4, 
p_0 : =y"5- (9/4) *y~4- (30/17) *y~3*x+(27/16) *y"3+(45/17) *y"2*x 

+(225/289) *y*x~2- (27/64) *y"2- (135/136) *y*x- (675/1156) *x~2 

£>5 being the common denominator, a conductor element lying in P, though 
AW = x 9 is the one computed directly from the Jacobian. The implementation 
also produces a Grobner basis B for the presentation: 

p_0~2- (135/17) *p_0+(81/4)*p_l*p_5~3-27*p_2*p_5~5-81*p_2*p_5~2 

-243*p_3*p_5"5- (405/17) *p_4*p_5"4- (243/8) *p_4*p_5"3 

- (1215/17) *p_4*p_5+ (729/4) *p_5"5 , 
p_0*p_l-9*p_0*p_5~2- (135/17) *p_l- (27/2) *p_2*p_5 

- (81/4) *p_3*p_5~4-27*p_4*p_5"6- (405/17) *p_5~5 

+(243/16)*p_5~4, 
p_0*p_2-27*p_l*p_5"4-81*p_l*p_5- (135/17) *p_2+ (81/2) *p_4*p_5"4 

+ (243/4) *p_4*p_5-243*p_5~6 , 
p_0*p_3-9*p_l*p_5- (15/17) *p_2+ (27/4) *p_4*p_5-27*p_5~6 , 
p_0*p_4-9*p_2*p_5-27*p_3*p_5"4- (135/17) *p_4+ (81/4) *p_5"4, 
p_l"2- (9/4) *p_0*p_5-9*p_l*p_5"2- (15/17) *p_2*p_5- (9/4) *p_2 

+ (27/4) *p_4*p_5~2-27*p_5~7 , 
p_l*p_2- (27/2) *p_l-9*p_2*p_5"2-27*p_3*p_5"5- (135/17) *p_4*p_5 
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+(81/8) *p_4- (81/4) *p_5"5 , 
p_l*p_3- (3/2) *p_l-p_2*p_5"2- (15/17) *p_4*p_5+ (9/8) *p_4 , 
p_l*p_4-p_0*p_5- (3/2) *p_2 , 
p_2"2-9*p_0*p_5- (27/4) *p_2-27*p_4*p_5~5 , 

p_2*p_3-p_0*p_5- (3/4) *p_2 , p_2*p_4-9*p_l+ (27/4) *p_4-27*p_5~5 , 
p_3~2- (3/4) *p_3-p_4*p_5~2- (15/17) *p_5 , 
p_3*p_4-p_l+(3/4)*p_4, 
p_4~2-p_2 

with induced weights wt(p ) = 25, wt{p\) = 21, wt{p 2 ) = 20, wt(p 3 ) = 11, 
wt{pi) = 10, and wt(p^) = 6. [Note that (po,Pi,P2,P3,Pi,P5) here corre- 
spond to (y 4 , y 3 , y 2 , y~i, x\) in the notation of this paper; but the notation 
C/*25) /2I) /20) /lO; /e) would be better than either, given that the subcripts 
then correspond to the weights.] 

Using Macaulay2's integralClosure function, [7], an implementation of de 
Jong's algorithm, [1], the output ideal is generated by: 

314432y6-8489664xll-2829888y3x4-707472y5-832320y4x+2122416y2x4 
+530604y4+2496960yx5+1248480y3x+734400y2x2-132651y3-468180y2x 
-550800yx2-216000x3, 

w_ (3 , 0) x2-68y2+51y+60x , 

4624w_ (3 , 0)y4-6936w_ (3 , 0)y3-8160w_ (3 , 0)y2x+2601w_ (3 , 0) y2 

+6120w_(3,0)yx-8489664x9-2829888y3x2+2122416y2x2 

+2496960yx3+244800y2-183600y-216000x, 
68w_ (3 , 0) "2y2-51w_ (3 , 0) "2y-60w_ (3 , 0) "2x-8489664x7-2829888y3 

+2122416y2+2496960yx w_ (3 , 0) "3-41616w_ (3 , 0) y-8489664x5 , 
w_ (4 , 0) x-4w_ (3 , 0) "2y+3w_ (3 , 0) "2 , 
17w_ (4 , 0) y-60w_ (3,0) ~2-41616w_ (3 , 0) yx-8489664x6 , 
w_ (4 , 0) w_ (3 , 0) -2448w_ (3 , 0) "2x- 146880w_ (3,0) -33958656yx4 

+25468992x4, 

w_ (4 , 0) "2- 146880w_ (4,0) -407503872w_ (3 , 0) y2+305627904w_ (3 , 0) y 
-9236754432y4x-83130789888yx5+20782697472y3x+8150077440y2x2 
+62348092416x5-15587023104y2x-12225116160yx2 
+3896755776yx+4584418560x2 I 

This is a presentation relative to Z, as attested to by the leading coefficients. 
One can use u^o = 686 and u^o = 17 * 68 2 a to clean this up a bit, but it will 
still be an affine 5*-algebra presentation with a default block ordering, grevlex 
on the new variables, forced to have the input over Z as an explicit subring. 
The fact that it essentially found a single common denominator conductor 
element in P is uncharacteristic. 

Trying Singular's normal function, [9] [3], also an implementation of de 
Jong's algorithm gives numerators: 
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68y2x3-51yx3-60x4 , 

4624y4x-6936y3x-8160y2x2+2601y2x+6120yx2+3600x3, 
18496y5-41616y4-32640y3x+31212y3+48960y2x+14400yx2-7803y2 
-18360yx-10800x2, 

x5 

The relations are: 

s [ 1] =314432*y~6-8489664*x~ll-2829888*y~3*x~4-707472*y-5 

-832320*y~4*x+2122416*y~2*x~4+530604*y~4+2496960*y*x~5 

+1248480*y~3*x+734400*y~2*x~2-132651*y~3-468180*y~2*x 

-550800*y*x~2-216000*x"3 

s[ 2]=68*T(3)*y~2*x-51*T(3)*y*x-60*T(3)*x~2-33958656*y*x~7 
-11319552*y~4+25468992*x~7+16979328*y~3+9987840*y"2*x 
-6367248*y"2-7490880*y*x 

s[ 3]=3468*T(3)*y*x"3-340*T(3)*y~2+255*T(3)*y+300*T(3)*x 
-1731891456*x~9+169793280*y*x~6-577297152*y~3*x~2 
-127344960*x~6+432972864*y~2*x~2+509379840*y*x~3 
+49939200*y~2-37454400*y-44064000*x 

s[ 4]=T(3)*x~5-18496*y~5+41616*y~4+32640*y~3*x-31212*y~3 

-48960*y~2*x-14400*y*x"2+7803*y~2+18360*y*x+10800*x~2 

s[ 5]=272*T(3)*y~3-2448*T(3)*x"4-408*T(3)*y"2-240*T(3)*y*x 

+153*T(3)*y+180*T(3)*x-135834624*y~2*x~6+203751936*y*x~6 

-76406976*x~6-39951360*y~3+59927040*y~2+35251200*y*x 

-22472640*y-26438400*x 
s[ 6]=60*T(2)*x-17*T(3)*y*x+8489664*x~7+2829888*y~3 

-2122416*y~2-2496960*y*x 
s[ 7]=4*T(2)*y-3*T(2)-T(3)*x 
s[ 8]=T(l)*x~2-68*y~2+51*y+60*x 

s[ 9]=41616*T(l)*y*x+60*T(2)-17*T(3)*y+8489664*x~6 

s [10] =4080*T ( 1) * y ~2-3060*T ( 1) *y-3600*T ( 1 ) *x-17*T (3) *y*x~2 

+8489664*x~8+2829888*y~3*x-2122416*y~2*x-2496960*y*x"2 
s[ll]=5*T(3)"2-1797811200*T(l)*x-8489664*T(3)*y*x"2 

-734400*T(3)-46183772160*y~4*x+4239670284288*x~8 

-415653949440*y*x~5+1517136915456*y~3*x+40750387200*y~2*x~2 

+311740462080*x~5-1137852686592*y~2*x-1308087429120*y*x~2 

+19483778880*y*x+22922092800*x~2 
s[12] =T(2)*T(3)-41616*T(3)*y-2309188608*y~3*x~2 

+3463782912*y~2*x~2+2037519360*y*x~3-1298918592*y*x~2 

-1528139520*x~3 
s [13] =5*T (1) *T (3) -734400*T (1) -3468*T (3) *y*x+1731891456*x~7 
-169793280*y*x~4 

+577297152*y~3+127344960*x~4-432972864*y~2-509379840*y*x 
s[14] =T(2)"2-31212*T(2)-10404*T(3)*x-577297152*y~2*x"3 

+432972864*y*x~3+509379840*x~4 
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s[15]=T(l)*T(2)-41616*T(l)*y-8489664*x~5 
s[16]=T(l)~2-T(2) 



Again, this is written relative to Z and can be cleaned up a bit by using 
T(l) = 68 * c, T(2) = 68 2 * 6, T(3) = 17 * 68 2 * a This is at least a strict affine 
5-algebra presentation, but again suffers from being relative to S and having 
no hint of the induced monomial ordering. The fact that it found a conductor 
element in P is uncharacterisitc. 

In Magma, [8], the Normalisation function, a third implementation of de 
Jong's algorithm, gives a basis: 

[$ . 1-4-4913/3375*$ . 1*$ . 3~6+4913/1500*$ . 1*$ . 3"5-4913/2000*$ . 1*$ . 3~4 
+4913/8000*$ . 1*$ . 3"3+17/405*$ . 2~4-4913/10125*$ . 2~3*$ . 3~3 
+4913/13500*$ . 2~3*$ . 3~2+289/6075*$ . 2~2*$ . 3*$ . 4-17/45*$ . 2~2*$ . 3 
+4913/91125*$ . 3~2*$ . 4~2-578/675*$ . 3~2*$ . 4 , 

$ . 1~3*$ . 3-3/4*$ . l"3-289/225*$ . 1*$ . 3"5+289/100*$ . 1*$ . 3~4 
-867/400*$ . 1*$ . 3"3+867/1600*$ . 1*$ . 3~2-289/675*$ . 2~3*$ . 3~2 
+289/900*$ . 2~3*$ . 3+17/405*$ . 2~2*$ . 4+289/6075*$ . 3*$ . 4~2 
-34/45*$. 3*$. 4 , 

$. 1~2*$. 2+15/17*$. l-$.3"2+3/4*$. 3, 

$. l~2*$.3"2-3/2*$. 1"2*$. 3+9/16*$. 1~2-17/15*$ . 1*$ . 3~4 
+51/20*$ . 1*$ . 3"3-153/80*$ . 1*$ . 3~2+153/320*$ . 1*$ . 3 
-17/45*$ . 2"3*$ . 3+17/60*$ . 2~3+17/405*$ . 4"2-l/3*$ . 4 , 

$. 1*$. 2"2+4913/1125*$. 1*$ . 3~7-4913/375*$ . 1*$.3"6 
+14739/1000*$ . 1*$ . 3"5-14739/2000*$ . 1*$ . 3-4 
+44217/32000*$ . 1*$ . 3-3-17/135*$ . 2"4*$ . 3+17/180*$ . 2~4 
+4913/3375*$ . 2"3*$ . 3~4-4913/2250*$ . 2~3*$ . 3~3 
+4913/6000*$ . 2"3*$ . 3"2-289/2025*$ . 2"2*$ . 3~2*$ . 4 
+17/15*$. 2~2*$. 3"2+289/2700*$. 2~2*$. 3*$. 4-17/20*$. 2~2*$. 3 
-1/9*$ . 2*$ . 4+15/17*$ . 2-4913/30375*$ . 3~3*$ . 4~2 
+578/225*$ . 3"3*$ . 4+4913/40500*$ . 3~2*$ . 4~2-289/150*$ . 3~2*$ . 4 , 

$ . 1*$ . 2*$ . 3~7-3*$ . 1*$ . 2*$ . 3~6+27/8*$ . 1*$ . 2*$ . 3~5 

-27/16*$ . 1*$ . 2*$ . 3"4+81/256*$ . 1*$ . 2*$ . 3~3-25/867*$ . 2~5*$ . 3 

+25/1156*$. 2"5+l/3*$.2"4*$.3"4-l/2*$.2"4*$.3"3 

+3/16*$. 2"4*$.3"2-5/153*$.2"3*$.3"2*$. 4+75/289*$. 2~3*$.3"2 

+5/204*$ . 2~3*$ . 3*$ . 4-225/1156*$ . 2~3*$ . 3-125/4913*$ . 2"2*$ . 4 

-l/27*$.2*$.3~3*$.4~2+10/17*$.2*$.3-3*$.4+l/36*$.2*$.3~2*$.4~2 

-15/34*$ . 2*$ . 3"2*$ . 4+1125/4913*$ . 3*$ . 4 , 

$ . 1*$ . 3" 14-6*$ . 1*$ . 3~13+63/4*$ . 1*$ . 3~12-189/8*$ . 1*$ . 3" 11 

+2835/128*$ . 1*$ . 3~10-1701/128*$ . 1*$ . 3~9+5103/1024*$ . 1*$ . 3"8 
-2187/2048*$ . 1*$ . 3"7+6561/65536*$ . 1*$ . 3~6 
+9375/1419857*$ . 2~6*$ . 3-28125/5679428*$ . 2~6 
-375/4913*$. 2"5*$. 3"4+1125/9826*$. 2"5*$. 3-3 
-3375/78608*$ . 2"5*$ . 3"2-25/867*$ . 2"4*$ . 3"8 
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+125/1156*$. 2~4*$.3~7-375/2312*$.2~4*$.3~6 
+1125/9248*$ . 2~4*$ . 3~5-3375/73984*$ . 2~4*$ . 3~4 
+2025/295936*$ . 2~4*$ . 3"3+625/83521*$ . 2~4*$ . 3"2*$ . 4 
-84375/1419857*$ . 2~4*$ . 3"2-1875/334084*$ . 2~4*$ . 3*$ . 4 
+253125/5679428*$ . 2~4*$ . 3+1/3*$ . 2"3*$ . 3~ll-3/2*$ . 2~3*$ . 3~10 
+45/16*$ . 2"3*$ . 3~9-45/16*$ . 2~3*$ . 3~8+405/256*$ . 2~3*$ . 3~7 
-243/512*$ . 2~3*$ . 3"6+243/4096*$ . 2~3*$ . 3~5 
+140625/24137569*$ . 2~3*$ . 4-5/153*$ . 2~2*$ . 3~9*$ . 4 
+75/289*$ . 2~2*$ . 3~9+25/204*$ . 2"2*$ . 3~8*$ . 4 
-1125/1156*$ . 2~2*$ . 3"8-25/136*$ . 2~2*$ . 3~7*$ . 4 
+3375/2312*$ . 2~2*$ . 3"7+75/544*$ . 2~2*$ . 3~6*$ . 4 
-10125/9248*$ . 2~2*$ . 3~6-225/4352*$ . 2~2*$ . 3~5*$ . 4 
+30375/73984*$ . 2~2*$ . 3"5+135/17408*$ . 2~2*$ . 3"4*$ . 4 
-18225/295936*$ . 2~2*$ . 3~4+125/14739*$ . 2~2*$ . 3~3*$ . 4"2 
-11250/83521*$ . 2~2*$ . 3"3*$ . 4-125/19652*$ . 2~2*$ . 3~2*$ . 4"2 
+16875/167042*$ . 2~2*$ . 3~2*$ . 4-125/4913*$ . 2*$ . 3~7*$ . 4 
+16875/83521*$ . 2*$ . 3~7+375/4913*$ . 2*$ . 3~6*$ . 4 
-50625/83521*$ . 2*$ . 3~6-3375/39304*$ . 2*$ . 3~5*$ . 4 
+455625/668168*$ . 2*$ . 3~5+3375/78608*$ . 2*$ . 3~4*$ . 4 
-455625/1336336*$ . 2*$ . 3~4-10125/1257728*$ . 2*$ . 3~3*$ . 4 
+1366875/21381376*$ . 2*$ . 3~3-1265625/24137569*$ . 2*$ . 3*$ . 4 
-1/27*$. 3~10*$.4~2+10/17*$.3~ 10*$. 4+5/36*$. 3~9*$.4~2 
-75/34*$ . 3~9*$ . 4-5/24*$ . 3~8*$ . 4~2+225/68*$ . 3~8*$ . 4 
+5/32*$ . 3~7*$ . 4~2-675/272*$ . 3~7*$ . 4-15/256*$ . 3"6*$ . 4"2 
+2025/2176*$ . 3"6*$ . 4+9/1024*$ . 3~5*$ . 4~2-1215/8704*$ . 3~5*$ . 4 , 
$.1*$.4 - $.2~2*$.3 + 3/4*$. 2~2] 

But the computed Grobner basis for this is way to big to be reproduced here, 
since Normalisation almost always chooses a default Lex monomial ordering 
the reverse of that which will give a readable Lex answer. As case-in-point, 
one can get a readable Lex monomial order answer by reversing the variables 
before computing a Grobner basis: 

$ . 1~2-135/17*$ . 1-$ . 3~5+9/4*$ . 3~3*$ . 4~3+27/4*$ . 3~3-405/17*$ . 3*$ . 4"4 

-243/16*$ . 3*$ . 4"3-243/4*$ . 4~8 , 
$. 1*$. 2-$. 3"3*$. 4-15/17*$. 3~2, 

$ . 1*$ . 3-27*$ . 2*$ . 4-4-9*$ . 3"2*$ . 4-135/17*$ . 3+81/4*$ . 4"4 , 
$ . 1*$ . 4-1/9*$ . 3-4+3/4*$ . 3-2+3*$ . 3*$ . 4-5 , 
$ . 2-2-3/4*$ . 2-$ . 3*$ . 4-2-15/17*$ .4, 
$.2*$. 3-1/9*$. 3-3+3*$. 4"5, 

$.2*$. 4-5-1/243*$. 3-5+1/36*$. 3-3+1/9*$. 3-2*$. 4-5+1/3*$. 3"2*$. 4-2 

+5/17*$ . 3*$ . 4-3/4*$ . 4-5 , 
$.3-6-27/4*$. 3-4-54*$. 3-3*$. 4-5-81*$. 3-3*$. 4-2-1215/17*$. 3~2*$. 4 

+729/4*$ . 3*$ . 4-5+729*$ . 4- 10 
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Since this implementation does not force a presentation relative to S, it oc- 
cassionally gives a decent minimized presentation. But there is, again, no hint 
that there is a natural induced monomial ordering. 



Since there is only one free variable in this example, Magma's IntegralClosure 
gives an answer 

B[l]= 1 
B[2]= Y 

B [3] = l/X~2*Y~2-3/4/X~2*Y-15/17/X 

B[4]= l/X~2*Y~3+(-15/17*X-9/16)/X~2*Y-45/68/X 

B[5]= l/X~4*Y~4-3/2/X~4*Y~3+(-30/17*X+9/16)/X~4*Y~2+45/34/X-3*Y 
+225/289/X~2 

B[6]= l/X"5*Y"5-9/4/X"5*Y"4+(-30/17*X+27/16)/X"5*Y~3 

+(45/17*X-27/64)/X~5*Y~2+(225/289*X-135/136)/X~4*Y 
-675/1156/X"3 

At least this necessarily gives a P-module basis and an answer over Q instead 
of Z. But there is obviously no way to give weights, and the presentation is 
only implicit. 
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